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CO  to 


ABSTRACT 


The  response  of  a  self-focusing  antenna  to  two  closely  spaced 
point  targets  is  investigated.  The  geometry  simulated  corresponds  to 
the  radar  tracking  of  a  low-flying  target  over  a  reflective  surface.  In 
one  case  the  antenna  consists  of  a  vertical  linear  array  of  elements. 
When  the  conjugates  of  the  incident  signals  are  retransmitted  at  each 
element,  after  several  iterations  the  antenna  radiation  pattern  intro¬ 
duces  significant  discrimination  against  the  weaker  of  the  two  unresolv- 
able  point  targets.  In  the  other  case  the  same  aperture  is  divided  into 
two  equal  subapertures,  and  conjugate  reflections  are  applied  to  the 
two  subapertures  rather  than  to  each  element.  It  is  found  that  very 
nearly  the  same  weak-target  discrimination  results.  The  implications 
of  these  results  with  regard  to  low-angel  tracking  are  discussed. 
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USE  OF  A  SELF-FOCUSING  ANTENNA  FOR 
LOW-ANGLE  TRACKING 


INTRODUCTION 

In  the  radar  tracking  of  an  object  positioned  over  a  reflective  surface,  a  reflection 
from  the  surface  may  appear  to  the  radar  to  be  an  image  below  the  surface.  (See  Fig.  1, 
taken  from  Ref.  1.)  If  the  energy  from  the  image  falls  within  the  main  lobe  of  the  track¬ 
ing  radar,  severe  errors  can  result.  There  have  been  a  number  of  attempts  to  reduce  these 
errors  (for  some  examples  see  Ref.  2),  but  the  problem  is  still  an  open  one.  This  report 
discusses  the  feasibility  of  using  an  adaptive  or  self-focusing  antenna  for  low-angle  tracking. 


Fig.  1  —  Reflection  from  the  surface  appears  as  a  target  image 

A  self-focusing  antenna  uses  the  unique  properties  of  conjugate  reflections  to  return, 
as  a  coherent  wavefront,  an  incident  wavefront  in  the  direction  from  whence  it  came.  If 
the  incident  waveform  is  backscatter  from  two  point  targets,  the  retransmitted  wave  is 
directed  toward  the  stronger  target,  thus  discriminating  against  the  weaker  target  (3-5). 

It  is  this  discrimination  against  a  weaker  target  which  could  render  this  technique  useful 
for  low-angle  tracking. 

At  the  time  this  investigation  was  begun,  it  was  well  known  (5)  that  discrimination 
against  a  weaker  target  would  occur  if  sufficient  aperture  were  available  to  resolve  the 
targets.  However,  if  self-focusing  is  to  be  useful  for  low-angle  tracking,  discrimination 
against  a  weaker  target,  i.e.,  the  image,  must  occur  when  the  target  and  image  lie  within 
the  main-lobe  beamwidth.  This  is  the  main  question  to  which  the  investigation  reported 
on  here  was  addressed. 
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MATHEMATICAL  DEVELOPMENT 
General  Case 


A  computer  simulation  program,  provided  as  an  appendix,  was  developed  to  calcu¬ 
late  and  plot  the  radiation  intensity  at  points  in  the  target/image  plane.  Plotting  the  ra¬ 
diation  pattern  proved  to  be  a  convenient  way  to  analyze  the  behavior  of  the  self-focusing 
antenna  for  closely  spaced  targets.  The  geometry  used  in  the  simulation  program  is 
shown  in  Fig.  2.  The  self-focusing  antenna  consists  of  a  linear  aperture  D  sampled  at 
every  half  wavelength.  The  sample  points  xi  are  a  measure  of  distance  above  the  refer¬ 
ence  plane  (ocean  surface).  The  number  of  aperture  samples  actually  used  in  the  program 
was  51.  The  distance  from  the  antenna  to  the  target  axis  is  R ,  and  the  target  height  is 
ht.  The  radiation  intensity  is  calculated  at  discrete  points  uf  along  the  target/image  axis. 
The  distance  dx y  from  a  sample  point  xi  along  the  aperture  to  a  sample  point  Uj  on  the 
target/image  axis  is 


d'j=s/R2  +(Xi-uj)  2.  (1) 

and  the  corresponding  electrical  length  of  that  path  is 

27 r 

tij  -y  dip  (2) 

where  X  is  the  wavelength.  Neglecting  all  attenuation  losses,  a  complex  transmission  co¬ 
efficient  r/;-  for  the  ijth  path  can  be  defined  as  the  complex  exponent  0{y,  i.e., 

r,y  =  j  >y,  (3a) 


where  N  is  the  total  number  of  sample  points  along  the  target/image  axis.  For  paths 
which  include  a  surface  reflection,  an  additional  180°  phase  shift  is  included,  i.e., 


r-  =  Jrtu  +  v  =  i  2 

1  ij  *  *  J  2  ’ 


(3b) 


The  attenuation  of  the  surface  reflection  could  be  accounted  for  in  the  coefficient 
however,  in  the  present  program  this  attenuation  is  accounted  for  by  reducing  the  cross 
section  of  the  image.  There  is  a  radar  cross  section  p(uj)  (a  complex  number)  associated 
with  each  sample  point  along  the  target/image  axis,  but  only  the  nonzero  points  corre¬ 
spond  to  the  target  and  its  image. 


The  process  is  begun  by  illuminating  each  point  along  the  target  axis  with  some  ini¬ 
tial  function  B0  (uj)  which  was  taken  to  be  uniform,  i.e.,  unity  for  all ;.  On  the  first 
iteration,  the  wave  back  at  the  ith  sample  point  in  the  aperture  is  given  by 

A  J  (Xi)  =  ^  B0  (Uj)  p  (Uj)  r,y.  (4) 

j 


The  paths  are  assumed  to  be  reciprocal,  that  is  Iy,  =  r,y.  At  each  point  xt  in  the  aperture, 
the  complex  conjugate  A*  (x,)  is  generated  and  transmitted.  This  is  the  essential  opera¬ 
tion  constituting  self-focusing.  In  practical  terms,  at  each  sample  point  in  the  aperture, 
the  phase  of  the  received  signal  is  reversed  and  then  retransmitted.  Back  at  the  target/image 
axis,  the  wave  at  each  point  is 
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Bx  (Uj) 


I 


A  I  (Xi)  Ty, 


(5) 


or  substituting  Eq.  (4)  into  Eq.  (5), 


B1  (“j>  =  ^ 


Bp  (Uk)  P*  ( uk )  F*k 

L  k 


(6) 


The  subscript  one  on  B1  (uj)  indicates  the  first  iteration.  After  iterating  this  process  a 
number  of  times,  results  converge  to  form  the  adapted  or  focused  pattern. 
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Fig.  2  —  Geometry  used  in  simulation  program 
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Case  of  Two  Subapertures 

A  variation  on  the  previous  approach,  which  can  lend  considerable  practicability  to 
the  concept  of  self-focusing,  involves  the  phasing  of  two  subapertures  rather  than  the 
phasing  of  each  element  of  an  array.  Two  subapertures  are  formed  by  dividing  the  aperture 
D  in  Fig.  2  into  halves.  Now  instead  of  retransmitting  the  conjugate  of  the  signals  received 
at  each  sample  point  in  the  aperture,  the  average  phase  across  the  aperture  is  computed, 
conjugated,  and  retransmitted.  On  retransmission  the  phase  across  the  subaperture  is  con¬ 
stant.  In  effect  then,  the  phase  distribution  across  the  entire  aperture  is  replaced  by  a 
two-step  approximation.  The  average  signal  for  the  first  subaperture  is 

M/2 

Ek  Ak  (Xi) ’  (7) 

i  =  1 
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where  M  is  the  total  number  of  sample  points  in  the  aperture  D  and  k  is  the  iteration 
number.  For  the  second  subaperture, 


M 


Fh  = 


Ak  (xi). 


■  M  , 

7  1 

Then  the  wave  at  a  point  on  the  target/image  axis  is  (feth  iteration) 

M/2  M 


Bk  E*k  r<7  F*k  ^ 


i=  1 


.  M 

,=ri 


M/2  M 

-«Sr"+F‘‘Sr«- 


‘“7  1 


(8) 


(9) 


Some  performance  degradation  is  expected  to  result  from  this  approach,  but  if  the  phase 
distribution  across  the  entire  aperture  is  benign  enough,  very  little  degradation  will  occur. 
This  is  the  second  major  question  to  which  this  brief  study  was  directed. 


RESULTS 

The  magnitude  J  Bk  (uj)  |  of  the  wave  is  computed  and  plotted  at  51  points  along  the 
target/image  axis.  (The  plot  routine  interpolates  linearly  between  computed  points.)  In 
Fig.  3  the  distance  between  the  target  and  its  image  corresponds  approximately  to  the  3-dB 
beamwidth  of  the  antenna.  The  arrows  at  1000  m  and  —1000  m  locate  the  target*  and  image 
respectively  and  also  indicate  their  relative  strengths.  A  distinct  increase  in  image  discrimina¬ 
tion  can  be  seen  for  the  second  iteration  (Fig.  3b),  and  for  the  fifth  iteration  (Fig.  3c)  the 
image  is  completely  nulled  out.  Positioning  of  a  null  exactly  on  the  image  is  not  possible 
in  general,  but  it  occurs  here  because  the  separation  between  the  target  and  image  corre¬ 
sponds  to  the  distance  between  the  peak  of  the  antenna  lobe  and  the  first  null.  The  pat¬ 
tern  has  converged  to  its  final  value  by  the  tenth  iteration. 

Figure  4  corresponds  to  a  half-beamwidth  separation  between  the  target  and  image. 

An  improvement  in  image  discrimination  can  be  seen  on  each  iteration  until  the  final 
value  is  reached,  A  close  examination  of  the  plots  reveals  that  little  main-lobe  distortion 
occurs;  however,  the  primary  effect  is  a  shifting  of  the  main-lobe  axis.  The  final  alignment 
of  the  radiation  pattern  is  such  that  the  target  is  strongly  favored  while  at  the  same  time 
the  image  is  illuminated  minimally. 

For  the  low-angle-tracking  problem,  a  more  realistic  ratio  of  target  strength  to  image 
strength  is  1  to  0.9.  This  case  is  plotted  in  Fig.  5  for  a  half-beamwidth  separation  between 
target  and  image.  Convergence  of  the  pattern  is  considerably  slower  than  in  Figs.  3  and  4 
where  the  relative  strengths  are  1  to  0.5.  By  iteration  30  the  discrimination  against  the 


♦For  an  S-band  (X  =  0.33  ft)  radar  and  a  6.6  ft  aperture,  this  corresponds  to  a  range  of  22  naut.  mi, 
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Fig.  3  —  Pattern  of  a  self-focusing  antenna  when  the  target/image  spacing  is  one  beam- 
width.  (a)  One  iteration,  (b)  Two  iterations,  (c)  Five  iterations,  (d)  Ten  iterations. 
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Fig.  4  —  Pattern  of  a  self-focusing  antenna  when  the  target/image  spacing  is  half  beam- 
width.  (a)  One  iteration,  (b)  Two  iterations,  (c)  Five  iterations,  (d)  Ten  iterations. 
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TARGET/IMAGE  AXIS  (m)  TARGET/IMAGE  AXIS  (m)  TARGET/IMAGE  AXIS  (m) 

Fig.  5  —  Pattern  of  a  self-focusing  antenna  when  the  target/image  spacing  is  half  beamwidth  and  the  relative  image  strength  is  0.9.  (a)  One 
iteration,  (b)  Two  iterations,  (c)  Five  iterations,  (d)  Ten  iterations,  (e)  Twenty  iterations,  (f)  Thirty  iterations. 
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O  -IOOO  0  500  1500  2500  -2000  -1000  0  500  1500  2500  -2000  -1000  0  500  I5C 

TARGET/IMAGE  AXIS  <m)  TARGET/IMAGE  AXIS  (m)  TARGET/IMAGE  AXIS  (m> 

Fig.  6  —  Pattern  for  self- focusing  applied  to  two  subapertures.  The  target/image  spacing  is  half-beam  width  and  the  relative  image  strength  is  0.9. 
(a)  One  iteration,  (b)  Two  iterations,  (c)  Five  iterations,  (d)  Ten  iterations,  (e)  Twenty  iterations,  (f)  Twenty-nine  iterations. 
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image  is  8.9  dB.  Although  not  shown  in  the  figure,  after  iteration  29  the  discrimination 
was  8.4  dB,  and  the  effect  is  somewhat  oscillatory.  The  final  value  would  be  about  8.7  dB. 
It  is  interesting  that  this  is  somewhat  better  than  the  final  value  in  Fig.  4d,  which  was 
6.75  dB. 

Finally,  in  Fig.  6  the  results  obtained  for  the  two  subapertures  are  shown.  The  re¬ 
sults  in  Fig.  6  are  extremely  encouraging  in  that  the  final  value  of  discrimination  is  about 
7.5  dB  or  only  about  1.25  dB  worse  than  the  results  for  the  full-aperture  case  of  Fig.  5. 
Convergence  is  slow  in  Fig.  6,  but  it  does  not  appear  to  be  any  slower  than  in  Fig.  5. 

There  is  some  cause  for  concern  in  Figs.  5  and  6  in  that  a  high  side  lobe  appears  below 
the  image  in  both  cases.  The  effect  of  this  side  lobe  is  to  increase  the  sea-clutter  level. 

It  is  possible  that  this  clutter  may  correspond  to  a  range  cell  outside  the  cell  of  interest, 
in  which  case  the  effect  may  be  negligible.  On  the  other  hand,  the  effect  of  the  clutter 
on  the  focused  pattern  itself  is  unknown. 


CONCLUSIONS 

It  has  been  demonstrated  that  a  self-focusing  antenna  can  discriminate  against  the 
weaker  of  two  targets  even  when  the  target  spacing  is  such  that  both  targets  are  well  within 
the  beamwidth  of  the  antenna.  Furthermore,  the  results  obtained  indicate  that  it  may  be 
possible  to  achieve  good  image  discrimination  in  low-angle  tracking  using  two  subapertures 
rather  than  phasing  the  individual  elements  of  an  array. 

The  results  of  this  brief  study  are  not  sufficient  to  conclude  that  a  self-focusing  an¬ 
tenna  would  be  useful  for  low-angle  tracking.  This  writer  feels  that  two  additional  ques¬ 
tions  must  be  considered  before  even  the  most  tentative  conclusion  can  be  drawn.  First, 
the  effect  of  target  motion  needs  to  be  considered.  This  can  be  investigated  by  a  rather 
straightforward  modification  to  the  existing  computer  program. 

The  more  critical  question  regards  the  actual  tracking  error.  To  investigate  this  ques¬ 
tion,  a  suitable  tracking  algorithm  must  be  developed  and  analyzed.  At  this  point  it  is 
not  perfectly  clear  how  to  use  the  self-focusing  technique  most  effectively.  However,  an 
approach  which  appears  promising  is  to  use  self- focusing  in  conjunction  with  monopulse: 
the  target  would  be  illuminated  by  the  self-focusing  antenna,  and  tracking  could  be  accom¬ 
plished  by  a  separate  monopulse  system  operating  in  a  receive-only  mode. 
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Appendix 


COMPUTER  SIMULATION  PROGRAM  FOR 
THE  SELF-FOCUSING  PROBLEM 


PROGRAM  C 

COMPLEX  C I ( 60 ) tCA ( 60 ) * A ( 60 ) tAO t GAMMA * DUMtZ ( 60 *60 ) *FAC<5*5) 

1  * AUX ( 60  *60 ) * RHO ( 60 ) *BUM 

DIMENSION  X ( 60 ) * U ( 60 ) *COSPH I ( 60 .60 ) tSlNPHI (60*60) *ANORM(60) • 

1AMAG ( 60 ) *NPLOT ( 20 ) *PLTARRAY ( 254 ) 

CALL  PL0T£(PLTARRAY.264*1) 

READ  2 *XMIN*XMAX*DELX 

2  FORMAT(6F10»2) 

READ  2  *UMIN .UMAX* DELL 
READ  3  *R  *NMAX  *  WLAM 

3  FORMAT (F10«2*I5*F10.2) 

READ  4 * ICASE  *NP * (NPLOT( I ) » 1*1 »NP) 

4  FORMAT (10(15)) 

ICASE°CAWE  NO.*NP°NO*  OF  PLOTS *NPLOT° I TERATON  FOR  PLOTS*NMAX«NO  OF  IT* 

Pl"3. 14159265 
AO° ( 1 • *0* ) 

GAMMA0 ( 1. *0* ) 

COMPUTE  X-ARRAY 

PRINT  5 • XMI N  *  XMAX  *  DELX 

5  FORMAT ( 1H1 . 15X.*XMlN°*.F10 .2 • laX **XM«X  =  *  *F 10 . 2  *  1 5X **DELX°* *F 10* 2 ) 

PRINT  6  *UMI N (UMAX  *DELU 

6  FORMAT ( //15X**UMIN°* *F10.2. 15X **UMAX°**F10.2  *  1 5X **DELU  =  * *F10 «2 > 
XLOD=WLAM/ ( XMAX-XMIN ) 

ELANG=UMAX/R 

XLUM=XLOD*R 

PRINT  7*R,NMAX*WLAM 

7  FORMAT (//15X**RANGE=**Fl0.2**Mc.TtRS  NMAX°* *  I  3 • 5X **LAMBDA=** 

1F6.2  »*CM* ) 

PRINT  8  *  XLOO  * LLANO *XLUM 

8  FORMAT  <//'10X.*LOD  =  *  .E  10.3  *5X**tL  ANGLE0* * E 10 .3  * 5X **LR/D=* * E 10 .3 ) 

I°1 

X(  I  ) °XM I N 

10  I-I+l 

X( I ) °X ( I“1 ) +DELX 
I F ( X ( I ) “XMAX ) 10*11*11 

11  IMAX=  I 

COMPUTE  U-ARRAYiJMAX 

NTARG=0 

J"1 

U ( J ) =UM I N 
RHO( J) °( 0*  *0*  ) 

12  J-J+l 

U(  J)=U(  J-D+DELU 
I F ( U ( J) .EQ.-500. ) 14.153 
153  IF ( U ( J ) «EO* 500* ) 15 *60 

15  RHO ( J )  =  ( 1 « *0.  ) 

JTARG=J 

GO  TO  16 

14  RHO ( J ) ° ( .9  *0  « ) 

J IMAG=J 

16  NTARG=NTARG+1 
GO  TO  17 

60  RHO ( J ) = ( 0* . 0 • ) 
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17  1F(U(J)“UMAX)12*13*18 
13  JMAX=J 

COMPUTt  S1NSCOS  OK  Pul 

JHALK=JKAX/2 
DO  19  1=1*1 MAX 
DO  19  J=1 *  JMAX 

Phl=Pl*2./WLAM*SURT(R**2+ (XI I >— U ( J1 1**2) 

IF ( J.LT.JhALF)PHl=PMl+Pl 
COSPH1  (  1  .  J  )  “COS  (  Prt  1  ) 

19  SlNPHl  (  i  *J  )  “Sli'KPnl  ) 

COMPUTt  COMP  tS  VALUtS  OF  UAMMA 

DO  20  1=1*1 MAX 
DO  2o  J  =  1  *  JMAX 

DUM  =  CMPLX  ( COSPFil  (  1  *J>  *S1NPII1  t  1  *J>  ) 

20  Z(  1  ,J ) =GAMMA*uUM 

CLEAR  AN ( U )  ARRAY  CNCt  ONLY  FOR  ALL  ITERATIONS 

DO  21  J  =  1  *  JMAX 

21  A(J)=(0.*0*) 

CLEAR  1N(U)  ARRAY  (TARGeT  REFLECTION)  AND  RECEPTION  ARRAYS  FOR  EACH 
ITERATION 

DO  22  J  =  1  * JMAX 

22  C1(J)=(0.*0.) 

DO  23  1  =  1  * IMAX 

23  CA ( 1 ) = ( 0 • • 0 • ) 

COMPUTt  FIRST  1TERAT ION  OF  AN ( U ) *  1 N ( U ) SAN ( X ) 

1 HALF 3 IMAX/ 2 
N=  1 

DO  24  J=1*JMAX 

24  C  1  ( J ) =RHO ( J ) *AO 
DO  25  1  =  1*1  MAX 
DO  25  J= 1  * JMAX 

25  CA(  1  >=CA(I  )+CK  J>*2(  I*J> 
tQ  1  IN  WATERS'  12 

BUM= (0.  *0* ) 

DUM= ( 0 • *0*  ) 

DO  128  1»1*1HALF 

128  DUM=DUM+CA ( I ) 

ME“ 1HALF+1 

DO  129  I»ME»IMAX 

129  BUM=BUM+CA ( I ) 

BUM=CONJG( BUM ) / ( IMAX-IHALF  ) 

DUM=C0NJG ( DUM) / 1HALF 

DO  130  1 “ 1  *  I  HALF 

130  CA (  1  )=DUM 

DO  131  1 “ME  * IMAX 

131  CA ( 1 )=BUM 

DO  26  J“1 * JMAX 
DO  26  1  =  1  * IMAX 

26  A(J)=A(J)+CA( 1 )*2( 1*J) 
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GO  TO  32 

CLEAR  IN(U)$AN(X)  ARRAYS  $  ITERATE 

37  00  27  J*  1 #  JMAX 

27  Cl ( J  )  s  (  0  •  #  0  • ) 

DO  28  1*1# IMAX 

28  CA ( I ) * ( 0 • *0 1 ) 

DO  29  J*1 # JMAX 

29  C I ( J)=KHO( J)*A( J) 

DO  30  1*1# IMAX 

DO  30  J  =  1 9 JMAX 

30  CA( I)*CA( IJfCI t J)*Z( I #J) 

BUM* (0- #0* ) 

DUM  = ( 0 • #0« ) 

DO  138  1*1# I HALF 

138  DUM=DUM+CA( I ) 

DO  139  I *ME  *  I MAX 

139  BUM* 6UM+CA ( I ) 

DUM*CON JG ( DUM ) / I  HALF 
BUM=CQNJG( BUM) / (  IMAX“ I  HALF  ) 

DO  140  I *1 # I  HALF 

140  CA ( I ) =DUM 

DO  141  I* ME# IMAX 

141  CA(  I  )=DUM 

DO  40  J - 1 » JMAX 
40  A( J)=(0.#0« ) 

DO  3 1 » J  =  1 # JMAX 
DO  31  1=1# IMAX 

31  A ( J ) =A ( J ) +C A ( I >*Z( I #J) 

32  PRINT  1003# I  CASE #NT ARG #N 

1003  FORMAT ( 14X#*ADAPTIVE  ARRAY  PHASING  PROBLEM*  # // # 17X • *CASE* # 1 2  #2X • 1 3 
1 #*P0INT  TARGERTS*»//20X#*ITLRATION  NUMBER* # I  3 # //13X #*U ( J ) * • 1 7X # 
2*NORMAL I  ZED  AN(U)*#//> 

AMAX  =  0 • 

DO  33  J  =  1 # JMAX 
AMAG ( J ) *C ABS ( A ( J ) ) 

33  AMAX=MAX1F ( AMAG( J) # AMAX ) 

NORMAL  I  Zt  AN  ( U )  *  PRINT  OUT  U  *  mN(U) 

DO  34  J* 1 # JMAX 

ANORM ( J ) *AMAG ( J ) / AMAX 

34  CONTINUE 

D I SCRM  =  2  0 • * A LOG 10  <  ANORM ( J  TAKu ) /ANORM ( J I  MAG ) ) 

PRINT  1003  #D ISCRM 
:.0J4  FORMAT  ( 11X  #F8.2# 10X  #l20»8 ) 

i.003  FORMAT  (  5X#*lMAot  JlSCRl.''ilNATiQi>i**#F6«2»*ED*) 

DO  50  I=1#NP 
IF(N.tu.NPLOT<  I  )  )  5 1 #  50 
51  CALL  GRAPHPLT ( ANORM # PLT ARRAY • JMAX ) 

50  CONTINUE 

IF(N-NMAX)35#36#36 

35  N*N  + 1 

GO  TO  37 

36  CONTlNUc 
CALL  STOPPLOT 
END 
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SUbROUT  INt  GRAPnPLT  (  aNOUki  *PLTArtrtAY  tJi-.AX  ) 
DIMENSION  AMKMC1)  t  PL  T  ARK  AY  ( 254 ) 

X  =  G  • 

Y*Q. 

IP— 3 

CALL  PLGT(X»Y»IP) 

I P  =  2 

DO  100  1*1*10 
X*X+  1 • 

CALL  PLOT ( X  t  Y  *  I P ) 

Y=  •  1 

CALL  PLOT ( X  #  Y  *  IP) 

Y*Q  • 

CALL  PLOT (  X  •  Y  » I  P  ) 

100  CONTINUE 
IP*  3 
Xs  0  • 

Y*G  • 

CALL  PLOT ( X • Y 1 1 P ) 

DO  200  1*1*10 

IP*2 

Y-Y+l. 

CALL  PLOT ( X  t  Y  *  IP) 

X*  •  1 

CALL  PLOT ( X  t  Y  » I P ) 

X*0* 

CALL  PLOT ( X  t  Y  t I P ) 

200  CONTINUE 

CALL  PLOT ( 4 • >0 • «  3 ) 

CALL  PLOT ( 4  • *9 • * 2  ) 

CALL  PLOT ( 4*  *0*  *2 ) 

CALL  PLOT ( 6 * #0.#3) 

CALL  PLOT (6*  *10* *2 ) 

CALL  PLOT ( 6*  #0*  *2  ) 

x*o* 

I P  *  3 

Y*  0  • 

CALL  PLOT  t  X  *  Y  t  I  P  ) 

DO  300  J= 1 1 UMAX 
I P  *2 

X*.2*( J-l) 

Y*10**ANORM( J) 

CALL  PLOT (  X  *  Y 1 1  P  ) 

300  CONTINUE 

CALL  PLOTS(0*0) 

CALL  PLOT S  (  0  *0  ) 

RETURN 

END 
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